Anti-Kondo regime of charge transport through a double dot molecule 
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The conductance through a serial double dot structure for which the inter-dot tunneling is stronger 
than the tunneling to the leads is studied using the numerical density matrix renormalization group 
method and analytic arguments. When the dots are occupied by 1 or 3 electrons the usual Kondo 
peak is obtained. For the case in which 2 electrons occupy the molecule a singlet is formed. Nev- 
ertheless, the conductance in that case has a constant non-zero value, and might even be equal to 
the maximum conductance of 2e 2 /h for certain values of the molecule parameters. We show that 
this is the result of the subtle interplay between the symmetric and anti-symmetric orbitals of the 
molecule caused by interactions and interference. 
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Double quantum dot devices have recently drawn much 
attention both for their relevance to possible technolog- 
ical applications such as qubits 0, 0, 0, 0], as we U 
as the light they shed on basic concepts such as the 
Kondo effect and the Fano resonance 

[Ell II 11 The 

double quantum dot molecule can be connected in par- 
allel to the leads, such that electrons may tunnel from 
each dot to any lead. In this case one expects inter- 
ference between paths going through each of the dots 
to play an important role in the transport through the 
molecule. The interplay between interference and inter- 
action effects such as Coulomb blockade, spin and or- 
bital Kondo effects, is the subject of many recent papers 
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In comparison, the transport through a serially con- 
nected molecule, for which electrons from the right lead 
can tunnel only to the right dot and vise versa, is ex- 
pected to be much duller. At very low temperatures one 
anticipates that for an odd filling of a molecule composed 
of two identical dots the Kondo effect will dominate the 
conductance. If the tunneling between the dots is much 
weaker than the coupling of the dot to the lead an or- 
bital Kondo effect (where the two degenerate states are 
an electron occupying the left or right dot) is expected. 
In the opposite limit a spin Kondo effect will be observed 
(where the degenerate states are a spin up or down elec- 
tron occupying the symmetric or anti-symmetric super- 
position of both dots orbitals). For a certain value of 
intermediate tunneling between the dots a SU(4) Kondo 
effect is expected. For even occupation of the molecule 
one may expect the usual Coulomb blockade scenario, 
i.e., once the number of electron on the molecule is inte- 
ger no conduction through it is expected. 

This naive picture may be broken by the following sce- 
nario: Each of the two electrons is localized on a sin- 
gle dot and forms a Kondo singlet with the electrons of 
the corresponding leads Transport is then possible 
between the left and right Kondo state, and the con- 
ductance may even reach the maximum value of 2e 2 /h 



for certain values of the molecule parameters [l£j . This 
state can be energetically favorable only if the gain from 
the formation of a Kondo state (which is of order of the 
Kondo temperature Tr- ) is larger than the loss of kinetic 
and interaction energy due to the localization of the elec- 
tron on a single dot (i.e., the triplet-singlet energy sepa- 
ration) . Since Tr depends exponentially on the coupling 
of the dots to the leads, one expects this behavior only if 
the hopping between the dot and the lead is much larger 
than the inter-dot hopping. 

In this paper we shall show that also in the opposite 
limit where the inter-dot hopping is much larger than 
the hopping to the leads, finite conductance for the dou- 
bly occupied molecule of two serial dots is possible and 
may even reach the maximum value of 2e 2 /h. Since in 
this state the ground state of the molecule is a singlet, 
this conductance is not connected to the Kondo effect. 
Rather, this conductance stems from interactions in the 
molecule, which result in partial occupation of both the 
symmetric and anti-symmetric states. Thus both states 
can carry current even when the average occupation of 
the molecule is integer. It should be noted however that 
this current vanishes in the limit when the symmetric and 
anti-symmetric states are degenerated - this limit corre- 
sponds to two disconnected dots. On the other hand in 
the opposite limit when the energy separation is large the 
situation is similar to the conventional Coulomb Block- 
ade, i.e., the current should be small for any integer occu- 
pation. All in all there should be some optimal splitting 
between the two states, which corresponds to a maxi- 
mum of the conductance. As we shall show the non-zero 
conductance valley is constant for a wide region of gate 
voltages and is only weakly dependent on the tempera- 
ture and external magnetic field, which is very different 
than the typical Kondo behavior. 

The serial double quantum dot molecule model is de- 
fined by the Hamiltonian 

H = ^molecule + H leads + -fornix- (1) 
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The Hamiltonian, -ff mo iccuic, of the double dot molecule 
is given by 

^molecule = ^2 (e + ha - ii)b\ a b ia + tb\ a b 2a + h.c 

i=l,2;<r=T,l 
i=l,2 

where the energy level e of each dot (with creation oper- 
ator &|(2)ct ^ or an e l ec t ron m the first (second) dot with 
spin a), the intra-dot hopping matrix element t and the 
charging energies U for each dot and the mutual capac- 
itive coupling energy U' . By applying a gate voltage V g 
to the molecule one may change the dot's level e to e — V g . 
In principle also the hopping between the dots, t, may be 
changed by applying a side gate (as well as the hopping 
between the dot and the lead). For any realistic configu- 
ration one expects U > U'. The influence of an external 
magnetic field h and chemical potential /x are included in 
the first term of Eq. J5J. The lead Hamiltonian, i?i ca ds — 
Ej>=T,J.,A=£,ii ~^]a\ a j^ + ( a } CT A a j+i<rA +h.c.) and the 
tunneling between lead and molecule is described by 
H mix = VY,a =hi {a\„L h ^ + a laR b ^) + t.c. Here at A 
is the creation operator in the i-th site of the A = L(R) 
left (right) lead, the hopping in the leads is set to one 
and the hopping between the lead and dot is V. 

It is useful to first consider the eigenvalues and eigen- 
vectors of a disconnected molecule, i.e., to diagonalize 
^molecule- Defining the single particle basis as |+) = 
(|1) + \2))/V2 and |-) = (|1) - |2»/V2 (where |1),|2) 
are the orbitals of the first and second dot), one can de- 
note the many particle Hilbcrt space by the application 
of the creation operators d) +a = (b\ a — b\ a )/\/2) and 
d}_ a = {b\ a — b\ a )l\p£) on the vacuum \vac). Thus, for ex- 
ample, | |, I) = d^d^^vac) and | tJ-jO) = d' + *(P + ^\vac) . 
The matrix corresponding to -ff molecule is block-diagonal 
with respect to the number of electrons ./V (./V may vary 
between zero and four) and each block may be diago- 
nalized independently. The following lowest eigenval- 
ues and eigenstates for each block are found: e(N = 
0) = 0,\vac) ; e{N = 1) = e - » - t;\ T,0) and 
e(N = 1) = e — fi — t; | J,, 0) (double degeneracy) ; e(N — 

2) = 2(e- M ) + (£/+£ /')/2 - y/(U - U'Y + (4t)72, (I + 

y i + [4t/(u - WW - 4t/(Z7 - uTr'il U,0) + 

Wl + \4tJJU - U')f - 4t/(U - U'))\0, U)} ; e(N = 

3) = 3(e - /j.) + U + 2U' - t; \ U, T) and e(N = 3) = 
3(e — /i) + U + 2U' — t; \ |J., J.) (double degeneracy) ; and 
e(N = 4) - 4(e - fi) + 2U + W, \ T|, T|>. 

Generally, one expects that coupling the molecule to 
the leads will cause broadening of the molecule states, 
but as long as V < t the state of the molecule will 
continue to retain their identity. Thus one might pre- 
dict the following behavior of the molecules occupation 
(i.e., number of electrons on the molecule) an function of 
e. For e > fi — t the dot population N = 0. Around 



e ~ /i — t the population switches to N — 1. At 
e ~ fx - (U+ U')/2 + ^(U-U') 2 + (4i) 2 /2 - 1 the popu- 
lation increases to N = 2, while at e ~ ^— U/2 — 3U'/2 — 
y/(U -U') 2 + (4t) 2 /2 + t it switches to N = 3. Finally 
at e ~ /i — U—2U' — t the molecule will be fully occupied. 
(2) According to the "orthodox" theory 13] one would 
a(T i = + \ expect a peak in the conductance each time when the 

population switches. Thus one would expect four peaks 
in the conductance as a function of e at the values 
of the switches in the molecule population given pre- 
viously. This consideration will not hold at zero tem- 
perature and zero magnetic field (T, h = 0), provided 
that Kondo physics will play an important role. Thus 
for odd occupancy of the molecule (N — 1,3), for 
which the ground state is degenerate, we expect to see 
a Kondo plateau in which the conductance is equal to 
2e 2 /h. On the other hand for the N = 2 occupancy 
there are no degeneracies in the molecule ground state, 
e.g., the singlet (S = 0) ground state is quite far away 
from the triplet state. Therefore naively one expects 
no conductance. Thus, at T = we expect two broad 
conduction peaks separated by a zero conductance val- 
ley. The width of the Kondo conductance peaks are 
'/m 2t+(U + U')/2- y/{U- U') 2 + (4i) 2 /2, while the width 
of the valley is -It + U' + yf{U -U') 2 + (it) 2 . 

Once we attach the molecule to the leads, the prob- 
lem becomes much more complicated. We compute the 
ground state for the interacting molecule attached to a 
couple of ID leads using an extension of a DMRG method 
previously used to calculate the ground state of a dot at- 
tached to leads JL4|. The essence of the method is sim- 
ilar to the regular DMRG for ID systems Once 
we obtain the many-body ground state eigenvector |0) 
of the entire system we can calculate the occupation 
of the |+) and |— ) orbitals of the molecule defined as 
n± a = (Q\d\ 1CT d± a \Q) . For the case of the double quantum 
dot with each of the two dots' is coupled symmetrically 
to a corresponding lead, and the interactions don't ex- 
plicitly break the symmetry, the Friedel sum rule implies 
the relation g = sin 2 (7r(n + | — n_f)) + sin 2 (7r(n + | — ti_j_)) 
in the usual way |l6j . Thus, the occupations n± a deter- 
mines the conductance. 

In Fig. we present the numerical results obtained 
for the conductance as function of the gate voltage (e) 
for different values of the inter-dot interaction U', while 
U = 0.8,t = 0.2 and V = 0.1 are kept constant. The 
two broad Kondo peaks (for which g = 2) are very clear 
and their width corresponds to our expectation 2t+ (U + 
U')/2 - ^{U - U') 2 + (4i) 2 /2 (see the inset of Fig. GJd) 
obtained for a disconnected molecule. The coupling with 
the leads rounds the peaks, but does not significantly 
change their width. Thus, in the absence of interactions 
(U, U' — 0) there is no peak (the width is equal to zero), 
while for U = U' the width is simply U. 

The surprise in Fig. ^ comes from the behavior of 
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FIG. 1: (a) Conductance, g, vs. gate voltage (level position), 
e, for different values of the inter-dot interaction U' , while 
the on-site interaction U = 0.8, inter-dot hopping t = 0.2 and 
dot-lead coupling V = 0.1 are kept constant and no external 
magnetic field is applied., (b) g, vs. e in the presence of a 
magnetic field h = 0.05 compared to g in the absence of such 
a field. Here U' = 0, and the other parameters are as in (a). 
Inset: the peak and valley width as function of [/'. Symbols 
correspond to the numerical results while lines to theory. 



narios. 

In fact a wide region of almost constant conductance 
for double dots connected in series could be seen in the 
data presented in studies using slave boson formalism 
|l7|. and numerical renormalization methods |l8|. Hints 
of this behavior are also present in the study for 3 dots 
in the N = 2 and N — 4 regions (lij . Nevertheless, 
the straight-forward expectation is that the conductance 
of the N = 2 region of the double dot molecule is zero 
unless some ferromagnetic coupling between the dots will 
create a triplet ground state for N = 2 [2(3, , or when 
each dot is more strongly coupled to its lead than to each 
other, in which case a separate left and right Kondo state 
form and transport between these states is possible 
Here none of these explanations is relevant. 

The mechanism behind the conductance plateau of the 
N = 2 state of the molecule has to do with the influence 
of the interaction on the ground state of the N = 2 sec- 
tor. In the absence of interaction, the + state (| Tl>0)) is 
full and thus can not conduct, while the — state (|0, tl)) 
is empty and thus also doesn't conduct. Due to interac- 
tions the molecule finds it favorable to occupy a super- 
position of both states and therefore neither of them are 
empty, and both, in principle may carry current. Since 
both states have an — 7r transmission phase difference, 
there will be interference between the two paths, which 
is taken into account in the Friedel sum conductance 
g = sin 2 (7r(n + | — n_|)) + sin 2 (7r(n + | — n_^)). Thus when 
71+ = n_ destructive interference will occur and the con- 
ductance will be zero, while when \n + — n_| = 1/2 maxi- 
mum conductance per spin of e 2 /h is obtained. Plugging 
in the values of n+ and n_ as function of t, U and U' 
previously calculated, we obtain the following expression 
for the valley plateau dimensionless conductance: 



the valley. Although its width corresponds quite well 
to the expectation of -2t + U' + ^{U - U') 2 + {At) 2 
(see the inset of Fig. [Qd) the conductance in the val- 
ley is constant and not necessarily equal to zero. The 
N = 2 state conductance is not suppressed by applying 
a magnetic field. As seen in Fig. \I]p, applying a mag- 
netic field that is enough to suppress the Kondo peaks 
of the N = 1 and N — 3 sectors, resulting in the ex- 
pected four conductance peaks, does not change much 
the conductance at the N — 2 valley. This is expected 
since the ground state of the disconnected molecule for 
the N = 2 sector is a singlet and does not couple to 
the external magnetic field. Only for a magnetic field 
h> (U'~ U)/2 + y/(U - U') 2 + (4t) 2 /2 will the triplet 
excited state (| ft, 0)) cross the singlet state and the con- 
ductance will drop to zero. Although the conductance in 
the valley has some superficial similarities to the Kondo 
conductance (i.e., is constant for a wide range of gate 
voltages), the fact that the ground state in the molecule 
has no degeneracies rules out any kind of Kondo- like sce- 
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7T(l-Wl + [4t/(U-U')] 2 -U/(U- 
1 + yi + [4t/{U - U')f - 44/(17 - 




This prediction was tested by comparing it to the nu- 
merical results for the conduction in the middle of the 
valley, depicted in the upper part of Fig. [2] A reason- 
able agreement can be seen. According to Eq. the 
conductance in the valley will be equal to its maximal 
value of 2 once the argument of the sin is equal to 7r/2, 
i.e., n+ — n_ = 1/2. Thus, for 



t 



U-U' 



(4) 



the valley conductance should be equal to 2 (if the 
molecule is asymmetricaly coupled to the leads, i.e., 
Vl 7^ Vr the valley height, as well as the Kondo peak 
height, will be 8{V L V R ) 2 /(V L 2 + V|) 2 ). As can be seen in 
the lower panel in Fig. for a double dot in which t is 
tuned to the value expressed in Eq. |@J, the conduction 
for all molecule populations ranging between N — 1 and 
N = 3 is equal to 2. The conductance value of g ~ 2 
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FIG. 2: (upper panel) Conductance at the middle of the valley 
for U = 0.8 and V = 0.05: as function of t while U' = (left) 
and as function of U' with t = 0.2 (right). The symbols 
corresponding to the numerical results fit quite well the line 
representing Eq. Q. (lower panel) Conductance as function 
of e for V = 0.1, U = 2, U' = 0.4 and t = 0.23 and no 
external magnetic field compared to the conductance in the 
presence of a magnetic field h = 0.05. With no magnetic field 
g ~ 2 for all molecule populations ranging between JV = 1 
and N = 3. When a weak magnetic field is applied gr ~ 2 
only for N — 2, while for a stronger field 4 Coulomb blockade 
peaks are observed. 



stems from two different processes. In the N = 1 and 
A^ = 3 regime it is the usual Kondo conductance while 
in the N = 2 it stems from the constructive interference 
between the conductance through the + and — states. 
This can be clearly seen when an external magnetic field 
is applied. The Kondo conductance is then suppressed, 
but the N = 2 state conductance isn't. 

Temperature is expected to have a very different in- 
fluence on the conductance for the different sectors of 
the gate voltage. As long as the molecule is in the 
Kondo regime (i.e., N — 1,3), the relevant energy scale 
is the Kondo temperature, Tk, which depends exponen- 
tially on the tunneling coupling of the molecule to the 
lead (V). On the other hand, in the "valley" (N = 2) 
as long as inelastic effects are ignored (in ID systems 
both electron-electron and electron-phonon scattering 
are rather weak) the conductance will change signifi- 
cantly only when ksT is of order of the singlet-triplet 



separation k B T ~ ([/' -U)/2+ ^{U - U') 2 + (4t) 2 /2. 

In conclusion, the conductance through a strongly 
bound identical double dot molecule was studied. The 
gate voltage applied on the molecule changes its fill- 
ing and therefore changes its conductance. Since the 
molecule is relatively weakly coupled to the leads, the 
number of electrons on the molecule is a good quantum 
number. For low temperatures (T < Tk) odd fillings 
(N = 1, 3) lead to the usual Kondo conductance of 2e 2 /h. 
For the even fillings of N = 0, 4 the conductance is equal 
to zero. The conductance for N = 2 is sensitive to the pa- 
rameters of the molecule, and can assume values between 
— 2e 2 /h. This behavior stems from the partial filling of 
both the symmetric and anti-symmetric orbitals due to 
interactions. Thus current fllows through both orbitals 
and the conductance is the result of the interference be- 
tween them. Thus, as function of the gate voltage a wide 
resonance tunneling plateau in the conductance appears, 
which is quite unsensative to temperature and applied 
magnetic field. 
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